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Abstract 


Stress  corrosion  cracking  and  general  corrosion  are  the  two  major 
processes  potentially  responsible  for  corrosion  failures  in 
prestressed  concrete.  Modeling  for  the  purpose  of  estimating  the 
corrosion  rates  of  steel  in  prestressed  concrete  are  discussed  in 
the  paper.  For  steel  in  concrete,  localized  corrosion  processes 
are  not  well  understood.  It  is  premature  at  present  to  develop  a 
sophisticated  but  still  incomplete  mathematical  model.  Instead, 
a simplified  approach  is  proposed  for  the  quantitative  evaluation 
of  corrosion  rates.  A model  was  developed  based  upon 
considerations  of  diffusion  processes  and  electrochemical 
reactions.  The  dominant  chemical  factors,  such  as  the 
concentration  of  oxygen  at  the  electrodes,  are  treated  as  primary 
variables  in  the  differential  equations. 

Two  approaches  which  have  been  frequently  used  to  treat  stress 
corrosion  cracking  are  briefly  described.  The  first  approach  is 
based  on  a conventional  engineering  approach  by  estimating  the 
stress  intensity  factors  under  various  environments.  In  the 
second  approach,  the  surface  energy  at  grain  boundaries  is 
estimated  along  the  crack  path.  It  is  suggested  that  the 
chemical  potentials  and  concentrations  of  the  important  absorbed 
species  be  included  in  the  energy  expression. 

Keywords:  concrete;  fracture  energy;  fracture  mechanics;  general 

corrosion;  modeling;  prestressing  steel;  stress  corrosion 
cracking . 
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1.  Introduction 


The  corrosion  of  steel  in  prestressed  concrete  (PC)  may  cause  a 
potential  threat  to  the-  integrity  of  the  structure  or  pipeline  of 
which  it  is  a part  even  though  the  concrete  in  PC  is  generally  of 
high  quality  (Treadway,  1971;  Szilard,  1969).  The  major  problems 
in  assessing  the  threat  have  been  identified  as  the  lack  of 
methods  for  estimating  the  corrosion  rates  and  for  evaluating  the 
effects  due  to  stress  corrosion  (e.g.,  Galligo  and  Climent 
(1980);  Wu  and  Clifton,  1981).  Because  of  the  complex  processes 
involved  in  the  corrosion  of  steel  in  concrete,  an  ad  hoc 
approach  was  proposed  in  reference  no.  4 for  computing  the  rates 
of  general  corrosion  processes.  Ordinary  differential  equations, 
with  time  as  the  independent  variable,  were  suggested  to  be  used 
in  making  the  computations.  This  paper  extends  the  previous 
effort.  A few  simple  parametric  functions  are  assumed  for 
illustration  and  solutions  are  plotted  and  discussed. 

For  stress  corrosion  cracking,  the  most  direct  treatment  of 
stress  corrosion  failures  appears  to  be  the  estimation  of  the 
fracture  stresses  in  terms  of  the  environmental  and  geometric 
parameters.  Semi-empirical  formulas  based  on  limited 
experimental  tests  are  traditionally  used  by  engineers.  A brief 
description  of  this  engineering  approach  is  given  in  this  paper. 
The  work-energy  concept  proposed  by  Cottrell  (1964)  appears  to 
provide  a more  fundamental  approach  that  can  correlate  failure 
mechanisms  with  materials  properties,  i.e.,  the  microstructural 
properties.  This  type  of  microstructure  approach  is  also 
described  in  the  paper.  It  is  suggested  that  the  expression  for 
the  surface  energy  include  the  chemical  potentials  of  the 
important  absorbed  species  at  the  grain  boundaries  of  the 
prestressing  steel. 
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2. 


General  Corrosion 


2 . 1 The  Analytical  Model 

Diffusion  and  electrochemical  reactions  are  the  major  processes 
involved  in  the  corrosion  of  steel.  Due  to  the  complicated 
conditions  existing  near  the  electrodes,  it  was  proposed  (Wu  and 
Clifton,  1981)  to  use  the  following  simplified  equations  in 
evaluating  the  rates  for  general  corrosion: 


^ - n ~ 
dt  “ ,2 


q (A*,  q)  q(qi) 


dVf 

dt 


= -Df 


+ q (A<|),  q)  + q(q,) 


(1) 


(2) 


Where  U^,  U^,  and  represents  the  concentrations  of  oxygen, 
ferrous  and  ferric  ions,  and  chloride  ions,  respectively,  near 
the  electrodes,  and  A4)  is  a term  representing  the  potential 
difference  between  electrodes.  The  terms  U*  and  denote  the 
concentrations  of  oxygen  and  of  ferrous  and  ferric  ions, 
respectively,  at  the  boundaries  of  ionic  and  hydrate  layers 
surrounding  the  prestressing  steel.  D^  is  the  equivalent 
diffusivity  of  oxygen.  Df  represents  the  diffusivity  of  ferrous 
and  ferric  ions  near  corrosion  sites  on  the  prestressing  steel. 
The  parameters  Aj  and  A2  link  the  thickness  of  effective  ionic 
and  hydrate  layers  around  the  corrosion  sites.  The  variables  Cj, 
C2,  C3  and  C4  are  functions  to  be  determined  from  the 
experimental  data.  In  an  exact  analytical  representation,  the 
ferrous  and  ferric  ions  should  be  considered  separately.  In  the 
approximate  representation  given  by  equation  (2)  their 
concentrations  are  combined. 

The  forms  of  the  above  equations  are  the  same.  They  may  be 
expressed  as: 


dt 


* 

aU^  = aU^  + b 


(3) 


with  Uj  representing  the  variables  or  Uf  in  equations  (1)  and 
(2) , and  a and  b representing  the  rest  of  the  variables. 
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2.2  Illustrative  Examples 

A few  simple  parameters  functions  are  assumed  for  a and  b as 
listed  in  table  1.  The  general  solution  for  equation  (3)  for 
these  cases  is: 


(b)  + K^ib) 


where  kj(b)  and  k2(b)  are  dependent  on  the  boundary  conditions. 
The  corresponding  solutions  of  kj  and  k2  for  the  assumed 
functions  are  also  listed  in  table  1.  To  illustrate  the  shapes 
of  the  transient  states,  solutions  of  U-  for  case  b in  table  1 
are  plotted  in  figure  1.  The  curves  b-1,  b-2 , and  b-3  represent 
solutions  for  U-  when  the  value  of  variable  a is  chosen  to  be 
0.1,  0.2,  and  0.5,  respectively . 

2.3  Discussion 

The  shapes  of  the  transient  states  from  solving  equation  (3)  are 
compared  with  those  from  the  traditional  approach  used  in  the 
field  of  electrochemistry,  i.e.,  characterization  of  the 
reactions  with  a set  of  partial  differential  equations  and 
solving  the  equations  using  conventional  techniques  such  as  the 
Laplace  Transform  method.  If  the  reaction  is  of  the  simple 
charge  transfer  type;  the  corrosion  current  can  be  expressed  as 
(e.g.,  Fraunhofer  and  Banks,  1972): 

at)  = s ^ (5) 

'ft 


where  S is  a coefficient  that  represents  the  effect  from  the 
significant  parameters  and  is  independent  of  time.  The  form  of 
i(t)/S  is  shown  as  curve  e in  figure  1 to  give  a comparison  with 
the  solutions  obtained  from  equation  (3) . The  limiting  states 
for  all  these  solutions  approach  zero. 
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3 . Stress  Corrosion  Cracking 

Traditionally,  the  macroscopic  approach  to  evaluation  of  the 
fracture  has  been  emphasized.  Descriptions  of  the  use  of  this 
engineering  approach  to  evaluate  the  stress  corrosion  cracking  of 
prestressing  steel  are  given  in  the  following  section.  The 
correlation  between  the  fracture  stress  and  microscopic  material 
properties  is  likely  to  lead  to  a more  realistic  solution  to  the 
problem.  Cottrell's  (1964)  microscopic  model  based  on 
dislocations  theory  has  the  most  general  theoretical  background. 

A description  of  this  approach  is  given  in  a later  section.  In 
the  case  that  the  specific  mechanism  is  known,  e.g.,  due  to  film 
rupture  or  hydrogen  embrittlement  (Pugh,  1977) , the  prediction  of 
failure  by  a mathematical  modeling  approach  is  simplified. 

3 . 1 Engineering  Approach 

The  most  direct  way  of  determining  the  fracture  stress  is 
probably  by  representing  the  stress  intensity  factor  (Kj)  by  the 
initial  crack  geometries.  For  plates,  Kj  has  been  conveniently 
expressed  in  terms  of  the  ratio  of  the  crack  length  a^  to  the 
width  w of  the  specimen  (Chu  and  colleagues,  1980) : 


ifj  = F ( -fi  ) 
^ w 


(6) 


The  function  of  F depends  on  the  material  properties  as  well  as 
on  the  geometry.  In  practice,  F is  usually  represented  as  a 
polynominal  in  terms  of  the  selected  variable: 


N 

= A 


ii) 

w 


(7) 


Where  a^  is  the  crack  width, 

w is  the  thickness  of  the  plate, 

Cj , nj , are  constants,  and 

N is  the  total  number  of  the  terms  used  to  represent  the 
function  F. 
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Table  1.  Solutions  of  U-  from  Equation  3 Based  on  Simple  Forms  of  the  Variables  a and 
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N was  taken  as  4,  and  nj  as  1/2  1,  3,  4,  successfully,  in  the 
consideration  of  stress  corrosion  cracking  of  steels  in  H2S 
solutions  (Chu  and  colleagues,  1980) . 

In  environments  which  may  assist  crack  growth,  the  magnitude  of 
critical  stress  intensity  factor  is  undoubtedly  reduced.  The 
amount  of  the  reduction  in  terms  of  the  yield  strength  is  usually 
estimated  based  on  the  corroded  conditions.  With  a similar 
approach,  the  stress  intensity  factor  K for  a prestressing  tendon 
or  wire  may  be  expressed  by  considering  pitting  conditions.  For 
example,  the  fracture  stress  condition  of  a prestressing  rod  was 
expressed  as  (Galligo  and  Climent,  1980;  Asitz,  1976); 

= o ( [-  . 8554  — - 0 . 16  42]  4 + 2 . 9 4 8 — ^ 0 . 3 3 09  ) (8) 


where  Kj^  is  the  critical  stress  intensity  factor  o is  the 
stress,  D is  the  diameter  of  the  prestressing  rod,  and  a and  b 
define  the  pit  (see  figure  2). 

In  studies  of  the  kinetics  of  corrosion  cracking,  the  fracture 
time  (tp)  is  usually  found  to  be  a function  of  the  stress 
intensity  factor.  The  t^.^  curve  generally  has  a different  shape 
with  different  test  media.  Two  common  shapes  are  shown  in  figure 
3 (Charbonnier  and  Margot-Marette , 1980)  with  t^  plotted  on  a 
logarithmic  scale.  It  is  not  difficult  to  find  appropriate 
mathematical  expressions  to  represent  these  data.  A similar 
approach  may  be  used  for  the  corrosion  of  prestressing  steel.  As 
stated  earlier,  this  engineering  approach  could  be  the  most 
direct  way  of  finding  the  time  to  failure. 
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Typical  stress 


Figure  2.  Schematic  of  Stress  Corrosion  Crack. 
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3.2  Metallographic  Approach 

Stress  corrosion  cracking  (SCC)  may  follow  either  an 
intercrystalline  or  a transcrystalline  path  (Logan,  1966) . For 
either  type  of  failure,  the  fracture  energy,  w,  can  be  expressed 
by: 

w = (9) 

where 

w^:  energy  related  to  the  stress  field  of  the  dislocation 
Wgi  surface  energy  of  the  new  cracks  in  the  crystal 
w^:  elastic  energy  of  the  crack 

w^:  work  done  in  increasing  the  volume  on  opening  up  the 
crack 

If  the  crack  growth  occurs  by  coalescence  of  dislocations  in  a 
body-centered  cubic  crystal,  the  expressions  for  the  variables 
are  given  by  Cottrell  (1964)  as: 


^^"4  4J? 

^ 4n{l-v)  ^ d 


= 2yd 


d^ 

® 8G  ^ 


ana^  d 


(9A) 

(9B) 

(90 


(9D) 


where  G is  the  modules  of  rigidity 

Y is  surface  energy 

V is  poisson's  ratio 

o is  the  applied  stress 

R is  the  effective  radius  of  the  stress  field 
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d is  the  crack  length,  and 
a is  the  lattice  constant 

a 

n is  the  dislocation  constant 

For  the  problem  considered  here,  the  surface  energy  term  needs  to 
be  reformed.  The  chemical  potentials  of  the  important  species, 
such  as  hydrogen,  must  be  considered  at  the  grain  boundary.  The 
following  modified  Gibb's  equation  is  proposed: 

dx  +2  |i.  dc.  + dE^  =0  (10) 

where  \i.  represents  the  chemical  potentials  for  the  species 
considered,  and  c.  represents  the  surface  concentrations.  The 
variable  dE^  is  introduced  to  compensate  for  the  discrepancies 
between  the  theoretical  equation  compared  to  the  actual  SCC 
behaviors.  The  value  of  dE^  should  be  determined  directly  from 
the  experimental  data. 

The  feasibility  of  linking  equation  (11)  with  those  for  general 
corrosion  (eqs.  (1)  and  (2))  to  mathematically  describe  the 
complete  SCC  process  will  be  determined  in  future  studies. 


11 


4 . Summary 

An  analytical  model  for  evaluating  the  corrosion  rates  of  steel 
in  prestressed  concrete  is  proposed.  As  a feasibility  study,  the 
transient  states  of  the  solutions  are  plotted  against  the  results 
obtained  from  a different  approach  often  used  in  the  field  of 
electrochemistry.  Both  approaches  may  lead  to  limiting  states  as 
time  tends  to  infinity.  For  stress  corrosion  cracking,  either 
the  engineering  approach  or  the  metallographic  approach  may  be 
used  to  estimate  the  fracture  stress  of  the  prestressing  steel . 
For  the  latter  approach,  the  surface  energy  is  linked  with 
chemical  potentials  and  the  concentrations  of  the  important 
absorbed  species  at  the  grain  boundaries.  The  Gibbs  equation 
with  a compensating  term  should  serve  such  a purpose. 
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